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ABSTRACT 


The heat transfer mechanism and temperature distribution in laser welding applications have a great impact on 
the quality of the weld bead geometry, mechanical properties and the resultant microstructure characterizations 
of the welding process. In this study, the effects of pulsed laser welding parameters including the frequency and 
pulse width on the melt velocity field and temperature distribution in dissimilar laser welding of stainless steel 
420(S.S 420) and stainless steel 304 (S.S 304) was investigated. Acomprehensive comparison was conducted 
through the numerical simulation and artificial neural network (ANN). The results of numerical simulation 
indicated that buoyancy force and Marangoni stress are the most important factors in the formation of the flow of 
liquid metal. Also, increasing the pulse width from 8 to 12 ms due to increasing the pulse energy, the temperature 
inthecenterofthemelt poolincreased about 250 °C. This leads toincreasing the convective heat transfer inthe 
molten pool and heat affected zone (HAZ). The temperature difference at a distance of 1 mm from the beam 
center at both metals at a frequency of 15 and 20 Hzis bout 58 and 75 °C, respectively. Furthermore, reducing 
the frequency to 5 Hz, due to diminishment of thermal energy absorption time, has clearly decreased the weld 
penetration depth in the workpiece. According to the ANN results, increasing both pulse duration and frequency 
has the significant effect on increasing melting ratio from 0.4 to 0.8 compared to the other input parameters. The 
ANN results confirmed that under the same input conditions, because of the differences in thermal conductivity 
coefficient, absorption coefficient and melting point of the two pieces, S.S 304 has experienced higher tem- 
peratures about 10% more than S.S 420. Also, among the 13 back propagation learning algorithms, the Bayesian 
regularization algorithm had the best performance. Among the number of different neurons in the hidden layer, 
comparison was performed to prevent network overfitting. The maximum relative error ofnetwork output data 
and target data for S.S 304 and S.S 420 temperatures and melting ratio were 7.297, 10.16 and 11.33%, 
respectively. 


1. Introduction 


because the ratio of heat dissipation in the depth of the piece to the 
width in this type of welding is high. This is due to the high density ofthe 


Connecting and cutting different parts is very important in all ther- 
mal industries. In recent years, various methods for connecting parts 
have been proposed. One of the popular methods in recent years has 
been the use of laser beams to weld and cut various materials [1,2]. The 
use of this technology has been considered by scientists and craftsmen 
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beam and its concentration in the welding area. These unique features 
have led tothe widespread use ofthis method in various industries such 
as automotive, aerospace, marine, construction and laboratory sciences 
[3-5]. One of the uses of laser welding is to connect dissimilar alloys to 
each other [6]. The use of this method has been reported in previous 
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Fig. 1. Schematic configuration of the laser welding process. 


Fig. 2. Uniform grid system for the calculation domain. 


studies due to the different properties of the materials [7,8]. The fusion 
zone contains elements present in both base metals. Therefore, itis very 
important to study the properties and changes based on laser parame- 
ters. Various numerical methods have been studied to evaluate the 
quality of welding. Also, the use ofnumerical simulation methods leads 
to increased accuracy and quality in the welding process and reduces 
production costs [9-12]. Numerical simulation ofthe welding process of 
twometals stainless steel 304 (S.S 304) and stainless steel 420 (S.S 420) 
due to the unique properties of these elements can lead to an ideal weld. 
Investigation of stainless-steel alloys due to good strength and corrosion 
resistance is of great importance. For this reason, ifan ideal weld is 
achieved, it can be used in the manufacture of industrial and medical 
tools and equipment. 


Inastudy, the results of their smooth particle hydrodynamics (SPH) 
simulations of Lu et al. [13] showed that the waveform size in the 
interface area increases with increasing laser energy and along the 
welding direction. This finding was correlated to normal stress, hori- 
zontal welding speed, and jet speed. In a numerical study, Lin et al. [14] 
investigated the dimensions of the molten pool and the temperature 
distribution trend in the welding process of two metals, stainless steel 
and copper. Their results showed that the shape ofthe heat field and the 
molten pool was asymmetric and deviated towards the AISI 304 sheet. 
They also showed that laser power is the most important determining 
factorin the dimensions ofthe molten pool. Inanumerical study, Lietal. 
[15] investigated the microstructure and mechanical properties of two 
metals, stainless steel and brass. Their results showed that laser energy is 
a very effective parameter in amplitude, wavelength, and bonding 
waveform. It was also found that the weldability between the flyer and 
the substance plate ata flight distance of 0.2 mm is the optimal process 
mode. Inalaser welding process between brass and 88308, Yu et al. [16] 
examined the welding properties. Their results showed that due to the 
low melting temperature and higher thermal conductivity of brass, an 
asymmetric melting poolhas been formed, in which the meltingofbrass 
has a greater share. Dongetal. [17] investigated the effect of process 
parameters on melting ratio and temperature distribution in laser 
weldingoftwoalloys, brassandS.$308. Theirresults showed that using 
this numerical method can reduce error, cost and time. Ding et al. [18] 
investigated the effect of process parameters on the shape of the weld 
and the temperature distribution in laser welding of stainless steel and 
brass alloys. Their results showed that the melt pool is asymmetric, 
which is due to the lower melting temperature and high heat transfer 
rate of the brass alloy, which eventually increases the melt volume of 
this alloy. In a laser welding process between austenitic stainless steel 
and copper alloy, Nguyen et al. [19] investigated the temperature field 
and the microstructure of the fusion zone. Their results showed that 
increasing the welding speed not only reduces the temperature but also 
significantly reduces the cracking in the fusion area. It was also found 
that increasing the power ofthe laser shifts the cracks in the fusion zone 
to the stainless-steel base due to the higher temperatures for the steel 


Table 1 
Chemical composition for 304 and 420 stainless steel alloys [35,36]. 
Composition C% Cr% Mn% P% Si% Ni% Fe% 
Weight AISI 304 0.08 18.4 1.06 0.03 0.34 8.9 Balanced 
AISI 420 0.15 13 0.17 0.04 0.46 0.13 Balanced 


Fig. 3. Double ellipsoid and cylindrical heat source. 


and copper. Saha et al. [20] optimized the laser welding process pa- 
rameters duringtheS.S316weldingusingthecentralcomposite design 
numerical method. Their results confirmed that laser power, focal dis- 
tance and welding speed are important factors for the quality of the 
weld. Prabakaran et al. [21] optimized the process parameters in the 
AISI316 and AISI 1018 laser welding process using Taguchi-based GRA 
numerical simulation. Their results showed that the parameters of laser 
power, welding speed, and focal length are effective in determining the 
quality of welding. Using this method, they confirmed that the welding 
performed after heat treatment and after welding has the best tensile 
strength. Ina numerical study, Huang et al. [22] examined the welding 
quality in the laser welding process of aluminum and steel alloys. Their 
results confirmed that by increasing the fluid flow behind the keyhole, 
the greater depth-to-width ratio can make the keyhole more unstable 
and porous. Bhatt et al. [23] investigated the effect of process 


parameters on the quality of laser welding performed between S.S 316 
and brass. The outcomes showed that with increasing energy and laser 
power, the penetration depthineach metalincreases, whichishigherfor 
brass alloys. 

In astudy, Kumar et al. [24] investigated the effect of the angle of 
incident of the laser beam on the workpiece surface in an austenitic 
stainless-steel welding process. Their results showed that the geometry 
of the weld poolis dependent on the deviation in the angle ofimpact and 
the shape of the weld pool at a high-level change from a hemispherical 
or almost spherical shape toa teardrop shape by decreasing the angle of 
impact. Chen et al. [25] investigated the effect of laser welding pa- 
rameters on the welding quality of stainless steel and copper. Their re- 
sults showed that the laser parameters have a strong effect on the 
dimensions of the HAZ. It was also found that copper melting reduces 
jointtoughness. Lietal. [26] investigated the effect oflaser power on 
dissimilar laser welding of brass and stainless steel. Their results illus- 
trated that when the laser power was less than 1846 W, the obtained 
weld has less defects and has better microstructure and mechanical 
properties. In anotherstudy, Lietal. [27] examined the quality ofthe 
weld between the H62 brass and S.S 316L alloys. Their outcomes illus- 
trated that the welded joints performed better with the brass overlap 
configuration on the steel and no intermetallic composition was 
observed within the joint. In a laser welding process, Sasaki et al. [28] 
investigated the welding quality obtained from brass and stainless steel. 
Their results showed that by shifting the position of laser radiation to 
brass, agood quality weld can beachieved. Dongetal. [29] examined 
thequality oflaser welding performed on two metals, copper and brass. 
The results showed that it is possible to achieve a penetration weld of 
suitable quality with the help ofthis method. It was also found thata 
suitable quality weld can be achieved with the help of optimal process 
parameters. Galun et al. [30] investigated the welding quality of stain- 
less steeland brass alloys inalaserwelding process. The results showed 
that the weld strength obtained from brass-steel is higher than brass- 
brass. It was also found that a weld strong enough to withstand the 
deep tensile process was obtained by this method. Zhang et al. [31] 
investigated numerically and experimentally the residual stress fields at 
different bending angles and evaluated the effect of welding paths on 
residual stress distribution. The findings show that choosing the right 
welding path can reduce residual stress values between 20 and 40 MPa. 
Gengetal. [32] conducted the dissimilar laser welding of S.S 304 and S. 
S 420 with pulsed Nd: YAG laser and obtained the temperature distri- 
bution, microstructure and mechanical properties of the welded zone. 
Theeffects of pulse width and frequency on the fusion zone temperature 


Table 2 
Thermophysical properties of stainless steel 304 and 420 [46-50]. 
Nomenclature Symbol Material Value Unit 
Density ofsolidphase ps AISI304 ps = 8020—0.501(T —298) kgm 3 
AISI 420 7860 
Density of liquid pi AISI304 pı = 6900—0.8(T—1727) kgm 3 
phase AISI 420 7000 
Solidus temperature Ts AISI 304 1673 k 
AISI 420 1727 
Liquidus temperature Tı AISI 304 1727 k 
AISI 420 1783 
Ambient temperature Ta 298 k 
AISI 304 k = 10.33 + 15.4 x 10 ?T —7 x 1077 T?k = 355.93 —196.8 x 298 < T < 16331644< T < 1672T > 1793 
Thermal k 10 3Tk = 6.6 + 12.14 x 10 ?T Wm kt 
conductivity AISI420 k = 20 + 61.5x 10- “Tk = 133.4—594.9x 10-*Tk =6.5+116.8x T< 17271727 < T < 1783T > 1783 
10~*T 
Specific heat c c 2 0.443 + 2x 10 ^T —8 x 10 1072 298 < T < 1727 Jg tk 
AISI304 
AISI 420 c = 1.92—1.587 x 10-?T—20.569 1150 < T< 1173T> 1173 
Melting latent heat La AISI304 2.90 x 10° Tkg- i 
AISI420 3.04 x 10? 
Stefan-Boltzmann o 5.67 x 10-8 Wm ?Kk-^ 
constant 
Dynamic viscosity u 0.0007 kgm !s | 
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Fig. 4. The architecture of the ANN. 


Table 3 


Mean relative error value for different ANN training algorithms [58,59]. 


Algorithm Abbreviation Mean Relative Error % 
S.S 304 S.S 420 Melting 
Temperature Temperature Ratio 
Levenberg- lm 2.264 4.496 2.766 
Marquardt 
Bayesian br 1.917 3.077 2.268 
regularization 
One-step secant oss 4.107 3.865 6.044 
Scaled conjugate scg 3.686 5.173 4.665 
gradient 
BFGS quasi-Newton bfg 4.152 3.928 5.441 
Conjugate gradient cgp 4.046 5.135 4.797 
backpropagation 
with Polak-Ribiére 
updates 
Conjugate gradient cgb 4.563 3.599 2.369 
backpropagation 
with Powell-Beale 
restarts 
Conjugate gradient cgf 2.555 4.427 4.978 
backpropagation 
with Fletcher- 
Reeves updates 
Resilient rp 2.792 4.147 5.409 
backpropagation 
Gradient descent gdx 3.718 5.278 5.164 
with momentum 
and adaptive 
learning rate 
Gradient descent gd 11.501 12.166 9.541 
Gradient descent gdm 13.481 11.088 8.870 
with momentum 
Gradient descent gda 4.516 4.351 5.080 


with adaptive 
learning rate 


Table 4 
Mean relative error value to estimate the number of neurons in the hidden layer. 


Neuron Mean Relative Error % 
S.S 304 S.S 420 Melting Ratio 
Temperature Temperature 

3 3.12 4.63 9.13 

4 2.98 4.23 8.86 

5 2.82 3.92 10.10 

6 3.77 4.71 5.88 

7 4.46 5.45 6.33 

8 4.20 5.70 6.77 


variation were investigated. Their results showed that the laser pulse 
duration and frequency had induced almost similar changes for heating 
and cooling cycles. 

The metal's microstructure is altered by the quick heating and 


cooling process, which also changes the metal's mechanical and physical 
characteristics. The material is subjected to thermal expansion and 
tensionasaresultofthehigh temperature, which can induce distortion 
and cracking of the weld. In order to ensure the desired qualities and 


structure ofthe welded material, temperature control during laser 
weldingis essential [33,34]. The main aim of this study is to investigate 


the effect of pulsed laser welding parameters, specifically pulse width 


and frequency, on the temperature and velocity distribution and the 
dimensionsofthemolten poolusingnumericalsimulations basedonthe 
finite volume method. Furthermore, ANN are employed to develop a 
predictive model for the molten pool dimensions based on the laser 
welding parameters. For predicting the temperature at the center of 
molten pool and the behavior of melt flow, finite volume method would 
be more appropriate. On the other hand, ANN can predict the temper- 
ature and approximate dimensions of melt pool at lower time and cost 
which in turn will be suitable for utilizing in industry. By utilizing ANN 
and the finite volume method, this research aims to provide a compre- 
hensive understanding of the thermal behavior during laser welding and 
to develop a useful tool for optimizing the laser welding process. 


2. Numerical simulation 


The temperature gradient and liquid metal’s flow in the molten pool 
Parane r a gitn tpesualit anc mechanical property Toe 


transfer and fluid flow can increase the welding quality and reduce the 
costs. In order to numerical simulation of continuous laser welding, a 
finite volume transient numerical model was used. A numerical code 
was also used to define the thermal model and the thermophysical 
properties with temperature variation. Fig. 1 showed a schematicofthe 
laser welding process and the dimensions ofthe pieces. In order to 
investigate the various laser welding parameters and validate the re- 
sults, thetime history oftemperature, accordingto Fig. 1, was obtained 
aroundthecenterofthe beam. According to the Fig. 2, the number of 
uniform grid elements after the grid study was considered to be 558748. 
Table 1 shows the composition of the elements in stainless steel 304 and 
420. 


2.1. Laser welding process modeling 


In keyhole laser welding, cloud of plasma is created at the surface of 
the molten pooland prevents the penetration ofthe partialoflaser beam 
into the keyhole. This causes that a part of heat input absorbed at the 
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Fig. 5. The network training process for a) temperature of S.S 304, b) temperature of S.S 420 and c) melting ratio. 


Table 5 
Theused laser parameters to study theeffectofpulse width. 
Sample Welding Pulse Frequency Current Focal 
Speed Duration (Hz) (A) position 
(mm/s) (ms) (mm) 
1 3.1 8 15 130 10) 
2 3.1 12 15 130 10) 


surface of the workpiece, which represents the inverse absorption of 
Bremsstrahlung and the rest of energy is absorbed by the absorption of 
Fresnel on the keyhole wall [37]. Therefore, the Gaussian elliptic and 
cylindrical heat source was used respectively for modeling the thermal 
energy absorbed at the surface and depth of the workpiece. Fig. 3 shows 
the schematic of the used thermalmodel. The equations for the elliptical 
and cylindrical heat source are presented in the Eqs. (1) to (3) [38,39]. 
Allequations ofthe thermal model were defined as a source term in the 
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In the above equations, q,is the cylindrical volume heat flux, qr, qf 
are elliptical heat source in the front and back of the laser beam, 
respectively, pis the laser power, f,, fare elliptical heat source distri- 
bution coefficients in the front and back ofthe laser beam, fi, f are the 
energy distribution coefficients for double ellipsoidal and cylindrical 
heat sources, respectively and a; a, b, c, d are heat source distribution 
parameters which is shown in Fig. 3. 


2.2. Governing equations 


In the simulation of laser welding, the convection heat transfer was 
occurred due to the liquid metalflow in the molten pool. Therefore, it is 


necessary to couple the momentum and energy equations with each 
other. To determine the melt and solid regions, the liquid fraction was 
defined as the following equation. When the meltis being freezing or 
solid is being melting, the amount of liquid fraction is between 0 and 1. 
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Fig. 6. Temperature history at pulse width of 8 ms at different distances from the center of the laser beam a) beam center, c) 1 mm, e) 2 mm, and pulse width of 12 
ms at different distances from the center of the laser beam, b) beam center, d) 1 mm, f) 2 mm. 
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Fig. 7. Temperature contour on the surface of the workpiece at pulse width of a) 8 ms and b) 12 ms. 


Table 6 


The used laser parameters to study the effect of frequency. 


Sample Welding Pulse Frequency Current Focal 
Speed Duration (Hz) (A) position 
(mm/s) (ms) (mm) 

3 6.2 9 15 130 0 

4 6.2 9 20 130 [0] 


Therefore, by adding a liquid fraction in the momentum equation, the 
force caused bythe change ofthe materialphase can be applied [38]. 
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The discretization of the equations was considered according to these 
assumptions [38-42]: 


- The flow of fluid was considered incompressible, Newtonian, and 
laminar. 

* The molten pool's surface was considered flat. 

* Vaporization of liquid metal was ignored. 

* The initial temperature of 25 °C was considered. 


Continuity equation: 
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On 


Where p, H, t, p, U,k, g, Kandy are density, total enthalpy, time, 


pressure, vector of velocity, coefficient of thermal conductivity, gravity 
acceleration, coefficient of Darcy resistance and dynamic viscosity, 
respectively [39-43]. 


2.3. Boundary conditions 


The boundary condition at the surface [44]: 
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Stefan-Boltzmann, Tis the surrounding temperatures, and his the 


coefficient of convection heat transfer. 
Shear stress due to surface tension gradient [45]: 
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The boundary condition on the lower and lateral surfaces was 
expressed d follow [44]: 
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2.4. Thermophysical properties 


Inlaser welding, the range of variation of parameters with temper- 
ature is high significantly. Therefore, some thermophysical properties in 
terms oftemperature were defined for pulsed laser welding simulation. 
Table 2 shows the thermophysical properties of S.S 304 and S.S 420. 


3. Back-propagation artificial neural network (BP-ANN) 


In this study, in order to foretell the effectoflaser pulse width and 
frequency on the temperature field and melting ratio, the ANN with the 
error propagation learning law was used. Neural networks with BP 
learning law are the most common types ofneural networks for function 
estimation and pattern recognition. In this type of network, a set of 
signals moves from input to output and other signals move from output 
to input, which are called Function Signals and Error Signals, respec- 
tively. In multilayer networks with BP learning law, the first category 
signals are calculated based on the function of each neuron input and 


network parameters, and the error signals change the network param- 
eters and adjust them to reduce the error [51-53]. 


k "T 7) A(T — To) + qs (8) 


3.1. Network architecture and training methods 


The feed network architecture intended for training consists of 
a hidden layer with a transfer function of sigmoid and an output layer 
with a transfer function of linear to convert the input signal to output 
signal. The presence of a hidden layer with functions which are 
nonlinear al- lows the network to apply well linear and nonlinear 


Where sis the coefficient of emission, o is the coefficient of 


relevance between input and output vectors. Therefore, the transfer 
function of tansig was applied in the hidden layer. This transfer 
function prevents the output values from increasing too much. The 
relationships of tangential and linear transfer functions are as 
follows [54]: 


2 
Yransig() = ; +1 (12) 
1-e7 


5500 5520 
5000 5000 
4500 4500 
4000 anao 
a520 Esso 
E a00 E smo 
= — 
£ 2500 § 2800 
2200 3 2000 
1500 1500 
1000 1900 
500 500 
a 4 35 5 $5 Us 4 a5 D $5 
Time (s) Time (s) 
c) d) 
700 700 
9600 500 
500 500 
g g 
g wo PTS 
s 2 
s = 
z xo & s 
E E 
200 200 
100 100 
qs : Ws 4 5 55 


Time {s} 


SS, M4 (Rep) 
5.5, 304 (Num,) 
SS, 420 (Exp) 
SS, 420 (Nera) 


45 
‘Time (5) 


a 5.5, M4 (Espo 
5.5. 9H (Num) 

(J 8.5. 420 (Expo 
— SS, (20 (Num) 


Temperature (°C) 


0 1 0 


10 
Time (5) 
Time (5) 


Fig. 8. Temperature history at frequency of 15 Hz at different distances from the center of the laser beam a) beam center, c) 1 mm, e) 2 mm, and frequency of 20 Hz 
at different distances from the center of the laser beam, b) beam center, d) 1 mm, f) 2 mm. 


Ypuretin(N) = n (13) biasisb, with the bias being provides atrainable constant value foreach 
neuron or node. 
In a neural network, the average relative error between the actual 


output and thenetwork outputs of the test datasets and trainingcan be 


Fig. 4 shows a view of the current network architecture whereyis the 
transfer function, the input vector is x, the weight vector is W and the 
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Fig. 9. Temperature contour on the surface of the workpiece at frequency of a) 15 Hz and b) 20 Hz. 
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Fig. 10. Velocity vector at the cross section of the workpiece at pulse width of a) 8 ms and b) 12 ms. 


a) 


Valocity 
? 220-01 


1719001 


$.5 420 


1 1450-01 


5731-02 


0.000e+00 
[mst 


b) 


2 8200-01 


2.122e-01 


7 07-02 


RA 


Fig. 11. Velocity vector at the cross section of the workpiece at frequency of a) 15 Hz and b) 20 Hz. 


Table 7 


The considered parameters for the validation ofthe molten pool [32]. 


Sample Welding Pulse Frequency Current Focal 
Speed Duration (Hz) (A) position 
(mm/s) (ms) (mm) 

8 3:3 13 18 138 =} 


used as acriterion to evaluate ofnetwork learning and performance 
[55]. Furthermore, the lower error value for both the training and test 
information is related to the better network performance. Also, mean 
square error (MSE) can also be used to evaluate the network perfor- 
mance [56,57]. In this study, the value of MSE in the error propagation 
algorithm was considered as the objective function and the network 
weights as thevariable function. The relationships between the average 


relative error and the average square error are as follows which n +y;and 


y; are the number of data, actual output and network output, respec- 
tively [58,59]. 


12 i— X 
Mean Relative Error = 7 ji Yel (14) 
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There are several algorithms for network training by Back- 
propagation method. 13 different algorithms [58-61] were considered 
in order to select the best algorithm for network training in terms of 
achievingthe leasterror. The average amount ofthe relative output 
errorforthevarioustrainingalgorithms was shownin Table3. As can be 
seen, the lowest error for all three parameters of temperature of S.S 304 
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Fig. 12. The molten pool's shape at welding speed of 3.1 m/s a) numerical simulation, b) real micrograph and welding speed of 4.3 m/s c) numerical simulation, d) 


real micrograph. 
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Fig. 13. The molten pool's shape at pulse width of a) 8 ms and b) 12 ms. 


and S.S 420 temperatures and melting ratio is 1.917, 3.0777 and 2.268 for 
the network with Bayesian regularization training algorithm, respec- 
tively. The Bayesian regularization algorithm does not require to a 
validation data setand the datais divided into only two categories of 
testing and training and this leads to the participation of all data in the 
training process [62]. 

This method mightstop the sooner convergences which reduces the 
error in each epoch. Therefore, Bayesian regularization learning algo- 
rithm was applied for network training. In determiningthe network 
architecture, the number of hidden layer’s neurons arevery important, 
because if the number of neurons is too small, the network will be 
underfit, andifitis toolarge, the dataused for training will be exactly 
overfit. But the network fluctuates sharply near these points. In this case, 
the network does not have the ability to generalize to the desired input. 
Therefore, to determine the appropriate amount of the hidden layer’s 
neurons, 3 to 8 neurons were considered. As can be seen in Table 4, the 
lowest errors were obtained in 5 and 6 neurons in the hidden layer for 
the temperature and melting ratio, respectively. 


Fig. 5 shows the performance of network training based on MSE, 
which was considered as the objective function to adjust weights and 
biases in the test and training phase in different epochs. According to 
Fig. 5, it can be seen that the lowest MSE in the training stage for tem- 
perature of S.S 304 and S.S 420 temperatures and melting ratio were 
50.87, 64.14 and 0.00078, respectively, and the network could not 
achieve lower MSE in subsequent epochs and training has stopped. Also, 
the network training for temperature of S.S 304, S.S 420 temperatures 
and melting ratio has stopped at 220, 1000 and 142 epochs, 
respectively. 


4. Results and discussions 
4.1. Numerical simulation 
Prediction of temperature distribution, velocity distribution and the 


dimensions of the molten pool in different input values can affect 
significantly the quality of weld. In general, according to different 
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Fig. 14. The molten pool’s shape at frequency of a) 15 Hz and b) 20 Hz. 
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Fig. 15. Network regression graph S.S 304 temperature for a) training data, b) test data and c) all data. 


thermalpropertiesofalloysindissimilarlaserwelding, determinationof 
the laser parameters range can be important to achieve the proper 
molten pool. For this purpose, in this section, the effect of pulsed laser 
welding parameters on the temperature and velocity distribution and 
the dimensions of the molten pool has been investigated. 


4.1.1. Temperature distribution 

Predicting the temperature distribution can help to control the pri- 
mary microstructure [63,64]. Hence, the temperature history was ob- 
tained at different pulse duration and frequency. 
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Fig. 16. Network regression graph of S.S 420 temperature for a) training data, b) test data and c) all data. 


4.1.1.1. Effectof pulsewidth. The pulse width has a significanteffecton 
the quality of weld, so that the created cracks caused by rapid freezing 
can be prevented with the correct adjustment of the pulse width [65]. 
Table 5 shows the parameters intended to investigate the effect of pulse 
width. 

In order to identify the effect of pulse width changes (8 and 12 ms) on 
temperature distribution, the time history of temperature at different 
transverse distances at 25 mm from the edge of the workpiece was 
investigated (see Fig. 6). According to the Fig. 6, it is observed that 
increasing the pulse width only increased the temperature at different 
times and had no effect on the temperature alternation. The temperature 
difference created between theS.S304.and8S.S420is different due tothe 
difference in their thermal conductivity coefficients. Due to the fact that 
the thermal conductivity coefficient of S.S 304 is less than S.S 420, the 
thermal penetration in S.S 304 is lower and the temperature near the 
laserbeamtobehigherthanofS.S 420. Butattimesthatthelaser beam 
has not yet reached to 25 mm from the edge ofthe workpiece or passed 
that,itisobservedthatthetemperatureofS.S 420is higherthanS.S304 
due to the low thermal penetration coefficient of S.S 420 at low tem- 
peratures. Fig.6.eand6.fshowacomparison ofthe temperature history 
obtained from numerical simulations and the experimental work which 
was performed by Geng et al. [32] at a lateral distance of 2 mm from the 
center ofthe laser beam. The results ofnumerical simulation arein good 
agreement with the experimental data which shows the correct defini- 
tion of thermophysical properties and thermal model. 

Fig. 7 shows the temperature contour in 15 mm from the edge of the 
workpiece at the end of the pulse width. It is observed that with 


increasing the pulse width due to the increasing the pulse energy, the 
temperature in the center of the melt pool increases by about 250 °C. 


4.1.1.2. Effect of frequency. The values intended to investigate the ef- 
fect offrequency on the temperature distribution, velocity distribution 
and molten pool dimensions are presented in Table 6. With the fre- 
quencyadjustmentinlaser welding, the heatinputtotheworkpiece and 
thermal cycle can be controlled. 

Fig. 8 shows the temperature history at frequencies of 15 and 20 Hz. 
It is observed that increasing the frequency caused to enhancing the 
temperature at various lateral distances. Also, increasing the frequency 
hasledtoincrease the temperature fluctuation at different times. As itis 
observed in the figure, the maximum molten pool temperature is ob- 
tained for both frequencies at 4.03 S, that with increasing the distance 
from the beam center due to heat penetration to the around, the 
maximum temperature is obtained at a higher time. The temperature 
difference between two pieces at a distance of 1 mm from the beam 
centeratafrequency of 15 and 20Hzisapproximately 58and 75 C, 
respectively. 

Temperature contour at surface of the workpiece in two frequencies 
of 15 and 20 Hzat 15 mm from the edge of the workpiece is shown in 
Fig. 9. Itis observed that due to the high thermal penetration coefficient, 
temperature distribution is diverted toward S.S 420. 


4.1.2. Velocity distribution 
Surface tension at the surface of the molten pool caused to the for- 
mation of the liquid metal flowin the molten pool. Created force caused 
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Fig. 17. Network regression graph of melting ratio for a) training data, b) test data and c) all data. 


by surface tension is changed with temperature and due to that the 
temperature is different at the molten pool surface, surface tension 
gradient is created at the surface ofthe workpiece and leads to the 
formation of Marangoni flow in the molten pool. 


4.1.2.1. Effect of pulse width. Due to the surface active elements in 
stainless steel, surface tension is reduced by increasing the temperature. 
Fig. 10 shows the velocity vector at the cross section ofthe workpiece at 
pulse widthsof8 and 12 ms. Itis observed that due tothe higher tem- 
perature gradientin S.S 304 comparing to S.S 420, the liquid metal flow 
is diverted towards the workpiece with a higher temperature gradient. 


4.1.2.2. Effect of frequency. The negative surface tension gradient 
caused to move the melt from the laser beam center to the molten pool’s 
wall. Increasing the shear stress at the surface in effect of the surface 
tension gradient leads to increase in the melt pool width. Fig. 11 shows 
vector of velocity at frequencies of 15 and 20 Hz. It is observed that with 
increasing frequency due to the increasing temperature gradient, liquid 
metal velocity increased. This leads to increasing the convective heat 
transfer in the molten pool and heat affected zone (HAZ). 


4.1.3. Shape and dimensions of the molten pool 

Predicting the dimensions ofthe weld bead in different laser pa- 
rameters can lead to the appropriate connection of the two sheets. In 
order to validate the results obtained from the finite volume method 
simulation of the molten pool, the micrographic results presented by 
Geng et al. [32] were utilized. Table 7 enumerates the parameters 
considered for validation. It is noteworthy that all the parameters are 


held constant, except for the welding speed which varies from 3.1 to 4.3 
mm/s. Notably, the focal position in the samples lies within 1 mm ofthe 
surface of the component. To account for the impact of the focal position 
on the simulation process, parameters m, c, and d, presented in equa- 
tions (1) to (3), were estimated to enable clear simulation ofthe effect of 
the focal position. 

Fig. 12 shows a comparison of the molten pool's dimensions and 
shape between actual micrographs results and numerical simulation. 
The figure shows how the simulated laser welding process successfully 
anticipated the size and shape of the molten pool, indicating that the 
thermal model's parameters were established. The outcomes showed 
that the numerical simulation captured the physical behavior of the laser 
welding process very well. As a result, the simulation approach can be 
regarded as a trustworthy tool for projecting how the molten pool will 
behave during laser welding. 


4.1.3.1. Effect of pulse width. With increasing the pulse width, the di- 
mensions ofthe melt pool and HAZ due to the rising the time of heat 
transfer were enhanced. Fig. 13 shows the liquid mass fraction at pulse 
widthsof8and 12 ms. Duetothelowthermal conductivity coefficient of 
S.S 304, the heat penetration to the around could be neglected. There- 
fore, the molten pool and HAZ dimensions for S.S 304 islarger than S.S 
420. 


4.1.3.2. Effect of frequency. Dimensions and formation of the molten 
pool have significant effect on the connection quality ofthe two pieces. 
As mentioned, due to different thermal properties of two workpieces, 
temperature asymmetric distribution was created at the surface and 
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Fig. 18. Relative Error diagram for a) S.S 304 temperature, b) S.S 420 temperature and c) melting ratio. 
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Fig. 21. Melting ratio estimated for different inputs: a) Welding Speed = 3.1 mm/s, b) Welding Speed = 6.2 mm/s. 


depth of the workpiece. Therefore, as shown in Fig. 14, the molten pool 
is asymmetric in the cross section of the workpiece. Itis observed that by 
reducing the frequency to the amount of 5 Hz, due to reduce of thermal 
energy absorption time, melt has not penetrated in the depth of the 
workpiece. 


4.2. Artificial neural network 


Figs. 15to 17 show the results obtained from the ANN in the training 
and testing phases. R correlation was determined to investigate the 
closeness of the output data from the neural network to the actual output 
in all diagrams. The line-passing equation of points was also determined. 
When two variables change in agreement with each other and in the 
same direction, their correlation coefficient is positive, and conversely, 
if two variables move in agreement with each other but in opposite di- 
rections, anegative correlationis created, andifthe variables do not 
agree to move with each other, there is no correlation between them. 
Correlation coefficient was used to measure the correlation between 
network output and actualoutput. For both datasets inthe trainingand 
test phase, thecloserthis coefficientis toone, thebetterthelearningand 


performance ofthe network. According to Fig. 15 upto 17, Rcorrelation 
fortemperatures ofS.S304,S.S420andmeltingratiowere0.964,0.944 
and 0.924, respectively. 

Fig. 18 shows the relative error ofthe actual output and the output of 
the trained network with respect tothe input parameters. It can be seen 
that the highest errors for temperatures of S.S 304, S.S 420 and melting 
ratio are 7.297, 10.16 and 11.33%, respectively which indicated the 
proper performance of the network in data training. 

As shown in Fig. 18, the network showed the best performance for 
the training data, but the network output is reliable when it also pro- 
vides logical outputs to the desired inputs. 

Figs. 19 to 21 show the output of the network in terms of desired 
inputs for temperature of S.S 304, S.S 420 temperature and melting 
ratio, respectively. The results indicated that the proposed model has 
satisfactory accuracy for allthree outputs. As shown in Figs. 19 and 20, 
under the same input conditions, because of the differences in thermal 
conductivity coefficient, absorption coefficient and melting point of the 
two pieces, S.S 304 has experienced higher temperatures than S.S 420. 
Also, according to the Fig. 21.c and 21.d, increasing the power of laser 
has the significant effect on the melting ratio compared to the other 


input parameters. 

Thejuxtaposition ofthe ANN outputs and the Finite Volume Method 
(FVM) simulation results in this particular research exhibited a 
commendable level of congruity between them. The predicted temper- 
ature history and temperature contours obtained from the FVM simu- 
lations were in good agreement with the ANN results. This serves as an 
indication of the robustness of the ANN in predicting the outcomes ofa 
system. In addition toits accuracy, theuse ofANN in comparison to FVM 
simulation can result in a noteworthy reduction in computational ex- 
penses. Furthermore, the ANN offers the capability to predict an 
extensive array of parameters, including the distribution of temperature 
as well as the dimensions ofthe molten pool. This makes it a highly 
versatileandefficienttoolformodelingandpredictingcomplexsystems. 
In comparison with other studies, the results obtained in this paper 
showed a similar trend in terms ofthe effect of pulse width and fre- 
quency on the temperature and the dimensions ofthe molten pool. 
However, it is important to note that the specific values of these pa- 
rameters and their effects may vary depending on the materials being 
welded and the specific experimental conditions. Therefore, it is 
necessary to carefully select the welding parameters and adjust them 
accordingly for each specific case to achieve the desired quality of the 
weld. A lot of attention has been given to dissimilar laser welding due to 
its potential for joining materials with different thermal properties. In 
this process, two or more materials with different compositions or 
structures are welded, which presents challenges in controlling the 
welding process and obtaining high-quality joints. 


5. Conclusions 


In this study, due to the expansion and need ofindustries to connect 
the dissimilar alloys, the effect of pulsed laser welding parameters on 
temperature, velocity and molten pool dimensions were investigated 
usingnumerical simulation and artificial neural network model. 

Generally, predicting the temperature at the center of molten pool 
and the behavior of melt flow and thermophysical properties which 
could be predicted effectively by finite volume method. ANN can predict 
the temperature mainly through the training process of data extracted 
form experiments and approximate dimensions of melt pool. 

In order to accurate modeling of laser welding, thermophysical 
properties were defined in terms of temperature and the elliptical and 
cylindrical heat sources were used. The results are as follows: 


. Due to the dissimilar laser welding of two workpieces, the temper- 
ature distribution was asymmetric and S.S 304 experienced a higher 
temperature than S.S 420 according to its lower thermal conductivity 
coefficient. 

: Increasing the pulse width from 8 to 12 ms, increased the molten 
pool’s temperature about 300°C. 

. Shear stress caused by surface tension gradient at the surface ofthe 
molten pool led to the flow of liquid metal in the molten pool. 

. With decreasing the frequency, the penetration of the melt in the 
workpiece was reduced. so that at frequency of 15 Hz, a proper melt 
pool was not created to connect the two pieces. 

. The presence ofa hidden layer with a sigmoidal transfer function led 
to the application of nonlinear relevance between output and input 
vectors to improve network performance. 

. The number of 5 and 6 neurons in the hidden layer created the best 
efficiency for network training for temperature and melting ratio, 
respectively. 

. The lowest MSEs in the training phase were obtained 50.87, 64.14 
and 0.00078 for temperature of S.S. 304, S.8 420 temperatures and 
melting ratio, respectively. 

. RCorrelation was obtained 0.964, 0.944 and 0.924 for temperature 
of S.S. 304, S.S 420 temperatures and melting ratio, respectively. 
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